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Abstract 



The charm quark offers interesting opportunities to cross-check the mechanism of CP violation pre- 
cisely tested in the strange and beauty sectors. In this paper, we exploit the angular and quantum 
correlations in the DD pairs produced through the decay of the V(3770) resonance in a charm 
factory to investigate CP-violation in two different ways. We build CP-violating observables in 
-0(3770) — > DD — > (V1V2XV3V4) to isolate specific New Physics effects in the charm sector. We also 
consider the case of ^(3770) — > DD — > (Vi^)(-K"7r) decays, which provide a new way to measure 
the strong phase difference 5 between Cabibbo-favoured and doubly-Cabibbo suppressed D decays 
required in the determination of the CKM angle 7. Neglecting the systematics, we give a first rough 
estimate of the sensitivities of these measurements at BES-III with an integrated luminosity of 20 
fb _1 at V , (3770) peak and at a future Super r-charm factory with a luminosity of 10 35 cm _2 s _1 . 
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1 Introduction 



Outstanding progress has been made over the last decade thanks to the data gathered at -B-factories, 
confirming that the Cabibbo-Kobayashi-Maskawa (CKM) mechanism embedded in the Standard 
Model (SM) is the main source of CP violation in the quark sector. The impressive agreement 
between results from the s-quark and the 6-quark sectors [U[2] calls for further checks in less tested 
areas. The recent discussions concerning the leptonic decays of D and D s mesons, about a possible 
disagreement between lattice results and experimental data [3H5], suggest that the charm sector has 
not been explored as extensively as other quarks [6j. Another illustration of this situation stems 
from .D-meson mixing, which has only very recently provided interesting tests of the SM and its 
extensions [7HTT]. 

Indeed, the D-meson sector is a remarkable place to improve our knowledge on CP violation in 
and beyond SM, for at least two different reasons. First, the SM predictions for CP violation in the 
charm sector are very small, due to the hierarchical structure of the CKM matrix and the difference 
of masses between the fermion generations. Any significant amount of CP violation would provide 
clear signals of New Physics, and a contrario, the absence of observation of CP violation already sets 
bounds on models beyond the SM [6j. Secondly, D decays play a prominent role in determining 7, 
the least well known of the three angles from the £>-meson unitarity triangle. A better understanding 
of the strong dynamics of related D decays would help in reducing the current uncertainty on this 
angle [EMU- 

On the experimental side, the final results from CLEO-c and the start of BES-III provide inter- 
esting opportunities. These charm factories are known to offer the possibility to exploit the quantum 
entanglement of DD pairs, as explained in several references [12 yi4ffT8] . In addition, it is also inter- 
esting to note that D — > VV (vector-vector) modes exhibit rather large branching ratios, of similar 
size with respect to the pseudoscalar (PP) or vector-pseudoscalar (VP) modes, and provide further 
angular observables to study the above issues. In this paper, we investigate this question which had 
not been detailed so far. 

In Section [21 we discuss production of coherent D°-D° pairs from ^(3770) decay, and in particular 
the angular distribution when at least one of the D meson decays into a pair of vector mesons. In 
Section [31 we apply these results to two different situations: the determination of CP-violating 
observables exploiting angular and quantum correlations in cases where both D decay into vector 
pairs, and the extraction of D — > Ktt hadronic parameters in relation with the measurement of the 
CKM angle 7. In Section HI we briefly discuss the application of these results for BES-III and Super 
r-charm factory, before concluding. 



2 Correlated D decays 

2.1 Basics 

We want to describe the decay chain for -0(3770) as 

V>(3770) -> D°D° -> (M 1 M 2 )(M 3 M 4 ). (1) 

where M\M 2 and M3M4 are mesons from two-body decays of D° and D°, respectively (Hereafter, 
ip denotes ■0(3770)). Since we do not tag the D mesons, and just observe their decay products, we 
can use two descriptions of D mesons - either the flavour states D° and D°, or the CP eigenstates 
(neglecting, for the sake of simplicity, CP-violation in D mixing): 

m Jj^m, mJ^m (2) 

with respective CP-parity: rjcp(Di) = —1 and r]cp(D2) = 1 (we take the convention CP\D°) = 
-\D°)). 
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Due to the spin of ip, the .D-pair is emitted with an orbital momentum L = 1 corresponding to 
an antisymmetric coherent state: 

|(DS)w) = d»». (3) 

One can in principle consider different situations as below, where V stands for a vector and P 
for a pseudoscalar meson: 

• (PP) + (PP), (PP) + (VP), (VP) + (VP): the only available observable is the branching ratio, 
since the partial waves and helicities are all fixed by angular momentum conservation. 

• (PP) + (VV), (VP) + (VV): (VV) can have three helicity states, and thus there are new angular 
observables. This can be exploited for (PP) = Kn in connection with the measurement of the 
CKM angle 7. 

• (VV) + (VV): this will be studied with an interest in new observables for CP- violation. 

The relevant modes for our studies can be extracted from Ref. [19] for the branching ratios and 
Ref. [20] for the projected efficiency at BES-III. 
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Table 1: Branching ratios for D decays into CP-eigenstates composed of two pseudoscalar mesons. In 
each case, the product of intrinsic CP parities and the estimated reconstruction efficiency at BES-III 
are indicated. Note that the efficiency is for both D° decaying into a PP final state; for single D 
decay, the efficiency is yfi. 
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Table 2: Branching ratios for D decays into CP-eigenstates composed of one pseudoscalar and one 
vector mesons. In each case, the product of intrinsic CP parities and the estimated reconstruction 
efficiency at BES-III are indicated. Note that the efficiency is for both D° decaying into a PV final 
state; for single D decay, the efficiency is y/e. 
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Table 3: Branching ratios for D decays into CP-eigenstates composed of two vector mesons. In each 
case, the product of intrinsic CP parities and the estimated reconstruction efficiency at BES-III are 
indicated. The rates in brackets are not measured yet, but were predicted in ref. [21]. Note that the 
efficiency is for both D° decaying into a VV final state; for single D decay, the efficiency is \fe. 



Some of the VV modes have not been measured yet, but some estimates combining naive fac- 
torisation and models for final-state interactions are available |21H24| . It is interesting to notice that 
the p + p~ decay mode has not been measured yet, even though one would expect it to be larger than 
p°p° (the latter being a colour suppressed mode). 

As it will become clear below, for the modes of definite CP-parity, one can sum over all the possible 
subsequent decay channels (enhancing therefore the branching ratios). For modes without such 
definite CP-parity, one needs to pick up a subsequent decay channel providing a definite CP parity. 
For instance, in the case of the modes where K*° is identified through the channel K*° — > Kstt , 
we have included a factor 1/6 in the branching ratio due to Clebsch-Gordan coefficients (l/y/3 from 
K*° to K°tt°, then l/y/2 from K° to K s ). 

Only the product of intrinsic parity is given, and one has to include the partial wave of the 
outgoing state. If PP as well as VV in S and D waves are not affected, one has to include an 
additional (—1) for P-wave states: 

Vcp(PV) = -r)cp(P)vcp(V), 7](VV,£ = 1) = -riCp(V)r,cp(V). (4) 



2.2 Differential decay width 

An adequate formalism to treat the question of decay chains is the framework of helicity amplitudes, 
described for instance in Refs. [25,26j. The decay chain is described by the product of amplitudes 
corresponding to each reaction. For a reaction A — > BC, we define polar angles (9a,<Pa) describing 
the momentum of particle B in the rest frame of A in a basis where the z-axis defined by the 
momentum of B in the rest frame of its mother particle. The decay amplitude depends on (9a, 4>a) 
and is denoted where A^, Ac* are the helicities of the daughter mesons. 

Let us start by describing the chain (we will study the other "path" later): 

i>^?D x D 2 , D x ^ViV 2 , D 2 ^V 3 V4, (5) 
Vi --> MiM{, V 2 --> M 2 M' 2 , V 3 -> M 3 M 3 , V 4 -»■ M4M4. (6) 

The helicity formalism yields an amplitude of the form 

A/rm _ A^D X D 2 a Vi-+M 1 M!l a V 3 ->M 3 M^ .V 2 ^M 2 M 2 V^AUM^ . Dl ^ Vl v 2 .d 2 ^v 3 v 4 , - 

m l2 — 00 ^00 ^00 ^00 ^00 A X Vl \v 2 Av" 3 \v 4 

Ay 

3E<o(^,^ 1 0)4 lDj (8) 
X V 
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xD% Xv3 _ Xv4 (ct> D2 ,9 D2 ,0)H% v ^ 



x 

r 3 9 




where m is the projection of the spin of the ip along an arbitrary axis, and Ay denotes collectively the 
helicities of the 4 vector mesons. The vector mesons are emitted from a spinless D-meson, so that: 
A = \ Vl = Ay 2 and k = Ay 3 = Ay 4 . We used the rotation matrix D J m , m (a, (3, 7) = e~ im ' a d j m , m {P)e~ imi 
with the Wigner d-matrix: 

d\ (9) = --LsmO, d 1 00 (9) = cos9, d 1 _ 10 (9) = ^=sm9, d° 00 (9) = 1. (9) 

The probability amplitude becomes 

MS = {^J^^^ 9 ^™^ ( 10 ) 

47r(47r) 3 e d m o(^)^ 

cos Vl cos 9 V2 Aq 1 ^' ViV2 - sin 0y, sin #y 2 cos <f> 12 A°^ VlV2 - sin 9 Vl sin #y 2 sin ^A^^ 2 

\/2 " \/2 

cos 0y 3 cos 9v 4 A° 2 ^ V3Vi - sin 0y 3 sin % cos S^A^ 2 ^ 14 - sin #y 3 sin 9 Vi sin $ 34 A^ 2 ^ V ' 3 14 

where we defined $12 = 4>Vi ~ 0V 2 an d ^34 = 4>v 3 ~~ 4>v 4 (i- e - the angle between the two relevant 
vector mesons) and the combination of amplitudes 

TT1IJV1V2V3V4 _ fJ"4> TtVi ttV 2 ttV A ttV 4 /-|-|\ 
n — D1D2 M\M' X M 2 M 2 M3M3 M4_M 4 ' 

A AA ' w V / 

We introduce the transversity amplitudes 

A|| = -^(ff+i + ff_i), Ao = i?o, A ± = -^(ff+i-fr_i). (13) 

My-, is actually only one of the two "paths" that can be chosen. The total amplitude for a given 
projection m of the spin of ip along an arbitrary z-axis is: M m = (-M$ + MR)/y/2. The differential 
decay width is obtained by averaging the squared modulus of the amplitude over the three possible 
values of m = +1,0,-1. The three squared Wigner functions (1^(9^) add up to 1, so that the 
differential width is 

81 

dT 4V = — — ,d(cos 9 Vl )d(cos 0y 2 )d$ i2 d(cos 9 V3 )d(cos 9 V4 )d^ 34 x \A^ v ^ Vi \ 2 (14) 

cos 9 Vl cos 9 V2 A°°^ VlV2 - — sin 9 Vl sin 9y 2 cos <f> 12 A°°^ VlV2 - -i= sin 9 Vl sin 9y 2 sin ^ 12 A^°^ VlV2 

V2 " v2 

cos 9 V3 cos 9v 4 A°°^ ViVi - — sin <9y 3 sin 9 Vi cos $uA^ V3V * - -^= sin 9 Vi sin 9 Vi sin $ 34 ^{ 14 

V2 " v2 

cos 9 Vl cos 9 V2 A°°^ VlV2 - sin 9 Vl sin % cos $ i2 ^f ^ VlV * - ~5= sin #y sin % sin $ i2 ylf °^ yiV2 

v2 v2 
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2 

cos V Va cos d VA Af^ VaV4 - sin 0y 3 sin 6 Vi cos $34^°^^ - -±= sin 9 Vi sin % sin $34^°^^ 

v2 " y2 

We have an integration over [0, w] for 0's and [0,2-7r] for and $34. The amplitudes A are 
normalized so that: T(X -> FZ) = -> FZ)| 2 . 

The above formalism can be adapted easily to describe the situation where one D meson decays 
into PP' rather than VV' . Indeed, it amounts to considering only the longitudinal decay amplitude 
for D — > V3V4 and to remove the angular phase space related to the decay products of V3 and V4. 

3 Observables from correlated D decays 

3.1 Observables from ip — > 2D — > AV for CP violation 

If we take the decay chain [15|.ll6j 

e+e- -> iji -> D°D° -> / a / 6 (15) 
with f a and /& CP eigenstates of same CP-parity, we have 

CP\^) = |V>, CP\f a f h ) = VaVbi-lYlfafb) = ~\faf b ) (16) 

since f a and f are in a P wave. Therefore, the decay of tp into states of identical CP parity is, by 
itself, a CP-violating observable pT5|[T6] . 

One obtains, neglecting CP-violation in DD mixing, the following result for the combined branch- 
ing ratio, which can be recovered from [T7] 

Bn((D D ) c= ^ -> f a f b ) = 2BK{D° -> f a )BTZ(D° -> / 6 )(|p a - Pb | 2 +r D |l - Pa/ 9 6 | 2 ), (17) 

with the ratio of CP-conjugate amplitudes and the combination of D-mixing parameters 

Pf = A(DoZ f f y rD = {x2 + y2)/2 < 10 " 4 - (18) 

where x = Am/T and y = Ar/(2r) are the difference of masses and widths of the mass eigenstates 
in the DD system, normalised by their average width [6]. 

If we assume that CP is conserved in decay, we have pf = rjf, and thus BIZ = for a, b with same 
CP-parity. Therefore, we have indeed that the observation of (D°D°)c=-i — > f a fb with a, b of same 
CP-parity is an indication of CP-violation. Let us notice that a and b must be different eigenstates 
(either different mesons, or for VV, different partial waves), and that this branching ratio is sensitive 
to different aspects of CP-violation compared to uncorrelated decays of D — > f a and D — > /&, since 
the latter would be sensitive to 1 — \p a \ 2 or 1 — \pb\ 2 - We can thus construct observables for CP 
violation in VV decays by considering states with the same CP parity, which depends on the relative 
angular momentum between the two mesons. The favourite channels among the measured ones are 
K+K-, 7T+7T-, K s tt ,p°it , K s p°, K*° p° -»• (K s tt°) (tt+vt^ ) and pV 

The transversity amplitudes A for D\ 2 — > W have simple transformation laws under CP: 

A n->vv> ^ + r, CP (V) VC p(V')ricp(D)A^ VV ' , (19) 
A^ VV ' -»• +Vcp(V)vcp(V')vcp(D)A^' , (20) 
A^ vv ' -> -Vcp(V)r,cp(V , )r jcP (D)A^ 9 ' . (21) 

Following eq. (]17p . CP conservation at the level of the amplitude would require that only two 
combinations of transversity amplitudes are allowed: (0, _L) or (||,_L). In terms of partial waves, 
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and || are combinations of S and D waves, whereas _L is P wave, which means that CP conservation 
at the level of the amplitude would impose the vector mesons to be emitted in (S, D) waves on one 
side and P wave on the other. 

Therefore, the following combinations of transversity amplitudes in the partial differential decay 
rate can be in principle CP violating observables: 

(0,0), (0,||), (11,0), (IUI), (j_,_L). (22) 

Let us notice that in the case of identical meson pairs in the final state, there is only one CP- 
violating configuration that is available: (0, ||), due to the Bose-Einstein statistics. It seems more 
interesting to consider two different meson pairs, both with longitudinal polarization, to get a larger 
BR. From the above table, the most interesting modes are K*°p° — > (i^57r°)(7r + 7r _ ) and p°p°. It is 
straightforward to construct the corresponding CP-violating observable: 



/ 



cffV— (5cos 2 6»vi - l)(5cos 2 6V 2 - l)(5cos 2 9 Vz - l)(5cos 2 #v4 - 1) (23) 
128 

_ I aWiV2V3Va\2\ 4-D°-*-ViV 2 |2| A D ^V 3 Vi\2 v i n „0 |2 

— I A I IA) I \ A o I x \Pvi,v 2 ~ Pv 3 ,vJ ■ 



Similar weights can be obtained for the other CP-violating combinations, exploiting orthogonality 
relationships for Legendre and Chebyshev polynomials to select specific angular dependences as in 
the previous case. For instance, we have: 



/ 



dr 4 y— (5cos 2 9 Vl -3)(5cos 2 6V 2 - 3) (5 cos 2 9 Va - 3) (5 cos 2 9 Vi -3) (24) 

,2 1 \( A „„„2 



•(4 cos $ i2 - l)(4cos $34 - 1) 
I aW\V2VzVa.\2\ a D*-^VJ/i\1\ a D°->V 3 V a \2 i l| || i2 

\ A I \ A \\ I \ A \\ I x \Pvi,v 2 Pv 3 ,vJ 



and 



/ 



dT 4V — (5 cos 2 (9^ -3)(5cos 2 6>y 2 - 3) (5 cos 2 9 Va - 3) (5 cos 2 9 Vi -3) (25) 
32 



•(4 cos 2 $12 - 3) (4 cos 2 $34 - 3) 

,v 2 ~ Pv 3 ,vJ 



|^ViV2F3F4|2|^D ^yil/ 2 |2|^D ^V 3 y4|2 x |^_L ,,1 2 



In the case of the same VV final state for both D decays, one can obtain the appropriate observable 
corresponding to (0, ||) using for instance 

J dT 4v[(5 cos 2 9 Vl - 1) (5 cos 2 6»y 2 - 1) (5 cos 2 9 Vs -3)(5cos 2 0vi - 3)(4cos 2 $i 2 - 1) (26) 
•[(5cos 2 f9 Vl - 3) (5 cos 2 6»y 2 -3) (5 cos 2 9 Vi - 1) (5 cos 2 9 Vi - l)(4cos 2 $ 34 - 1)] 

, _. _ „n T - w t _ T-^n , T T P _ II _ 



|^VlV 2 V3V 4 |2|^D ^V 1 V' 2 |2|^|D ^V' 3 y4|2 x |^0 



■ 2 

v 2 - Pv 3 ,vJ ■ 



3.2 yj ->■ 2D ->■ (VV)(Kir) for the extraction of 7 

The measurement of 7 from the Atwood-Dunietz-Soni (ADS) method |27| requires the determination 
of the hadronic parameters r and 5. At BES-III, we can also take advantage of the coherence of the 
D° mesons produced at the i/;(3770) peak to extract the strong phase difference S between doubly- 
Cabibbo-suppressed and Cabibbo-favoured decay amplitudes that appears in the 7 measurements [T2| 
Here we introduce, in the standard phase convention where 5 vanishes in the SU(3) limit, 

6 " (K-tt+\DV)- 10 
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The process of one D° decaying to K~ir + , while the other D° decaying to a VV CP eigenstate can 
be described as (for our purposes, it will prove more convenient to express this decay rate in terms 
of D° and D° amplitudes): 



dT 2V = ^-d(cos9 Vl )d(cos9 V2 )d<5> x\A^ VlV2 \ 2 \A D °^ Kw \ 2 
An 

cos9 Vl C os9 V2 (A?°-* VlV2 ~ re iS A^ v ^) 
~ S m9 Vl S m9y 2 cos - re iS A° 

-^sm9 Vl S m9 V2 sin^f ^ - re^f ^ 2 ) 
V2 



(28) 



We can introduce: 



a ,\\A d ° v *Vb) = MA D ° v a v b )p ' 11 



Va,V b - 



In the absence of CP violation, which we will assume in this section, we have: 

Pv!!,V b = -VGp(Va)vCp(V b ) = ~Pv a y b - 



(29) 



(30) 



Moreover, we notice that all the decays presented in Sec. 12.11 have CP parities such that p° = —1, 
which yields the further expression of the differential decay width in Eq. (128p 



dT 2V = ^-d{cos6 Vl )d{cose V2 )d<I>x\A^ VlV2 \ 2 \A D0 -> K7r \ 2 (31) 

An 

x cos 2 6 Vl cos 2 9 V2 \Af^ VlV2 \ 2 (l + 2r cos 5 + r 2 ) 

+- sin 2 9 Vl sin 2 9y 2 cos 2 $\A°°^ VlV2 \ 2 (l + 2r cos 5 + r 2 ) 

-V2 cos 9 Vl sin 9 Vl cos 9y 2 sin 9y 2 cos $Re[^°^ FlV2 (Af 0- ^ 2 )*] (1 + 2r cos 5 + r 2 ) 
+- sin 2 Vl sin 2 9 V2 sin 2 $| A^° 1 2 (1 - 2r cos 5 + r 2 ) 
+v / 2cos 0^ sin^x cos#y 2 sin0y 2 sin^jRelA^^Af^^ 2 )*]^ sin 5) 
+Im[<°^^(Af ^)*](1 - r 2 )} 

- sin 2 9 Vl sin 2 0y 2 cos $ sm${Re[A°°^ VlV2 (Af^ VlV2 )*}(2r sin 5) 
+Im[Af °^™(Af -^ 2 )*](1 - r 2 )}" . 

We see that the differential decay width provides six different angular observables depending on 
the following (real) quantities: 

• three products of moduli for VV decays: \A^ V *A D °^ K * A°^iV 2] 2 

• two relative phases between the three amplitudes Aj^,,^ 1/lV2 

• two strong parameters describing the Kn decay: r and 5 

Whereas the full angular integration yields the sum of the three transversity amplitudes: 



/ 
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dT 2V = ^\A^ V2 \ 2 \A D0 ^\ 2 \—\Af^ V2 \ 2 (l + 2rcos5 + r 2 ) + 
An L 9 



(32) 



9 1 II 



>ViV 2 \2 



(1 + 2r cos 5 + r 2 
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•Viy 2 i2 



(1 — 2r cos 5 + r 2 



one can easily separate the different contributions by choosing suitable weights for the angular integra- 
tion (they can be obtained easily by exploiting orthogonality relations among Legendre polynomials). 
In practice the best way to perform the experimental analysis is usually to do a maximum likelihood 
fit on Eq. (j3ll . 



Pi(e Vl ,e V2 ,^) 
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|(5cos^ -l)(5cos0f> 2 -1) 




M \ z 
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5^(5 cos Q\ - 3) (5 cos 9\ - 3) (4 cos $ 2 - 1) 




M„| 2 
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— cos dy 1 cos 6v 2 sin 9y 1 sin 9y 2 cos <3? 


Re 


M M* 





4 
5 
6 



"16 (J) COS V Vl 

^ cos 9 Vl sin 

1 /_ „o 



Vi cos 

4:V ^ _ J- _ 

■(5 cos g| - 3)(5cos6' 2 /ci - 



'y 2 sin #y 2 s i n 3* 
^ — 3) cos <I> sin 



|iW ± |" 

-Re [M M* 
Re [MjjMt] 



Table 4: Weights used to select contributions from the transversity amplitudes for ip — > (Ktt)(VV). 
M amplitudes are defined by Eq. ([33|h 



The branching ratio only depends on the three amplitude combinations: 



M = A, 



D°->ViV2 




(1 + re 



Mn = ^ 



d°^ViV 2 



(1 + re 1 



M, 



,4 



.D°->ViV 2 , 



1 



re 



(33) 



Table [5] shows that Pi. 2, 3 yield the relative size and phase of M$ and Mm, whereas Pi^fi yield the 
relative size and phase of Mj_. Therefore, without further knowledge, one can extract a combined 
constraint on r, 5 and A_\_ from the differential decay width, and one can also determine the ratio 
A|,/A . 

Note the invariance under the simultaneous transformation ^4o^.*| ~~ ^ ^o^lh ^4||^4j_ — — ^4*|^-_L> 
5 —7- —5, which implies that for fixed value of r and A± there is a twofold ambiguity on 5 (in other 
words there is no information on the sign of sin<5). It is worth noting here that the decay rate is 
sensitive to |sin<5| terms (thanks to Re [M 0i ||M^]), while in the standard analysis with PP modes 
one is only sensitive to cos 5 (neglecting the small mixing contributions which lift the ambiguity [12J ) . 
Since 5 is small, the sensitivity on the sine in addition to the cosine is expected to improve the final 
result. 

The above constraint can be improved by exploiting our current or expected knowledge of the 
polarisation of D — > VV . If we extract the relative size and phase of the three amplitudes A®^ ~^ VlV2 
from independent single D — >■ V1V2 decay (single-tag - ST) measurements, and if the three amplitudes 
are not too different in size (as seems to be the case for p°p°), the measurement of the Mj amplitudes 
in the correlated (double-tag - DT) DD — > (ViV2)(Kn) decay leads to the determination of both r and 
5 (more precisely, r, cos<5 and | sin <5| ). Since the ratio r is already well known, r = 0.055 ±0.002 [28], 
our method may lead to a good measurement of 5. 

Note that for relatively low statistics, a simplified transversity analysis can be performed. Instead 
of considering the full angular distribution in both single and double D decays, one can perform a 
one-parameter fit to the distribution of the transversity El angle 9t r , which yields the perpendicular 
polarisation fraction in single-tag and double-tag decays: 

' 2 |M ± | 2 



/ 



ST 



\A ± , 

A | 2 + |vl||| 2 + |^j 



/ 



DT 



|M | 2 + |M||| 2 + |M ± P 



9vi , 8ti, $tr) transversity angles are related to (9v 1 , &v 2 , helicity angles by 

cos 9v 2 = sin 9 tr cos 9 tI sin 9y 2 sin $ = cos 9 tI sin 9v 2 cos $ = sin 9 tI sin $ tr 



(35) 



(34) 



8 



The above observables leads to: 

\A ± \ 2 / ST 

which implies a one-dimensional parabolic constraint in the (r, cos 5) plane. An independent con- 
straint comes from the ratio of double-tag to single-tag widths proportional to (|Ao| 2 + |^4||| 2 + 
|-<4±| 2 )/(|-^o| 2 + l-^||| 2 + l-^±| 2 )> that can be expressed in terms of r, cos 5 and /^ T . This simplified 
transversity analysis allows one to determine r and cos 5. However the main novelty of our proposal 
comes from the sensitivity of the complete correlated decay rate to | sin <5| terms, which needs the 
study of the full angular dependence. 

3.3 CP-violation in D°D° mixing 

In the previous discussions, we have neglected the tiny CP-violation in D°D° mixing in order to 
simplify the study of correlated D — > VV decays. The inclusion of this effect would impact our 
results in the following way: 

• If CP violation is indeed measured through ^>(3770) — > D°D° — > (F1V2XV3V4), we cannot a 
priori disentangle CP violation in mixing from CP violation in decay. Therefore, if we want 
to convert this result into a bound on fundamental parameters of a New Physics model, we 
will have to exploit external inputs on CP-violating parameters of the mixing (from other 
observables). On the other hand, such an input is not necessary if we only aim at setting a 
constraint on CP-violation itself. 

• In the determination of the strong phase in D — > Kir, the amplitudes exhibit in principle a 
small dependence on mixing effects. However, this dependence is very weak with respect to the 
dependence on the hadronic parameters (r, 5), and as a first approximation, it can be neglected. 
As in the previous case, we can use external information on the CP- violating parameters in 
mixing to include their impact when required by more accurate measurements of the partial 
decay rate. 

4 Potential for BES-III and a super r-charm factory 

In this section we give a first rough estimate of the expected sensitivity of the two different measure- 
ments discussed above, either at the BES-III experiment or at a Super r-charm factory. 

4.1 CP violation 

As discussed in Section 13.11 the decay chain of e + e~ — > tj) — > D°D° — > f a fb can be described by 
Eq. (|17p . in which both CP conserving and violating processes can occur. We parameterize the ratio 
of amplitudes pj in Eq. (PT8|) as pj = r]j(l + 5f)e ia f , where the 5f is term from CP violation in decay, 
and aj is the phase difference between D° and D° decay into the same final state /. 

The D° decay channels in Table Q] can be directly used to search for CP violation by fully 
considering the correlation of D°D° production at BES-III. The background is small, and the main 
dilution is due to the mis-identification of charged particles, which is suppressed by about 10~ 4 . The 
sensitivity of measurement of CP violation can reach about 10 -3 with a 20 fb _1 luminosity on the 
^(3770) peak at BES-III. As in the previous case, one must take care of the background due to 
the dilution from non-CP eigenstates that impact the quasi two-body decays of D° meson listed in 
Table [2] (for D PV) and [3] (for D -> VV). 

Final states consisting of two vector meson pairs are particularly interesting, since one can use 
information on transversity amplitudes to extract different combinations of CP-violating observables, 



1 + re 



id 



1 



re 



is 



I"/ 



pDT 



(36) 
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as discussed in Section [3711 (0,0), (0, ||), (||,0), (||,||), (_L,-L). For example, a back-of-the-envelope 
computation yields the most promising channel p° p° / K*° p° : 

BK{{D°D°)c=-i P°p ,K* p ) CPV ~ 8 x BK\D° -> • -> ^*V) sin 2 ^=^, 

(37) 

where BIZ means the branching fraction for longitudinal polarized D° — > p°p° decay, BV} ' means the 
parallel helicities fraction of D° — > K*°p° decay, and where we have assumed that the CP-violating 
parameters 5f vanish. 

Assuming that no CP-violating signal events in D°D° coherent decays are observed with 20 
fb _1 data at BES-III, we can provide an upper limit on the CP-violating branching fraction at 90% 
confidence level (C.L.), as indicated in Table [5j A Super r-charm factory with 2 ab" 1 data yields 
naturally stronger constraints. If each polarized fraction is measured independently, an upper limit 
on the phase difference \a a — ctb\ can be set. For example, the current values for the polarized 
fractions in pp and pK* yields the upper limit \a a — aj| < 4.4° at 90% confidence level from the 

channel (D°D°)c=-i — > p°p ,K*°p ^ . At a future Super r-charm factory, with a data set of 2 

ab _1 , the constraint would be more severe, \a a — a&| < 0.5° at 90% confidence level. 

A more realistic analysis requires a likelihood fit to the full angular dependence of the VV modes. 
Systematics will arise from the mis-reconstruction as VV CP-eigenstates of the events that actually 
come from other resonances or background contributions. In view of the sizable width of the vector 
resonances, we expect that these systematics will dominate the final result. Their precise estimate 
in the framework of each experiment is however beyond the scope of this paper. 





Reaction 


Efficiency 


Upper limits 
at BES-III(xlO- 7 ) 


D»D" - 


-> (p+p-)(K*>u) 


0.13 


2.46 


D°D°- 


-> (p°p )(K*°p ) 


0.17 


1.88 


D°D°- 


+ {K* p )(K* uj) 


0.10 


3.19 


D°D° 


-»■ (K*°p ){p°<j)) 


0.09 


3.55 


D°D° 


-> (K^coKp ^) 


0.08 


3.99 


D°D° 


-»• (p°p°)(K*°uj) 


0.15 


2.13 


D°D° 


(pVXpV) 


0.13 


2.46 


D°D° 


^ (p + p-)(pV) 


0.11 


2.90 


D°D° —?■ 


(p+p-)(K*+K*-) 


0.11 


2.90 



Table 5: The projected 90%-C.L. upper limits on CP violating branching fraction of some most 
interesting (VV)(VV) modes from correlated D°D° pairs with 20 fb _1 data taken at ■0(3770) peak 
at BES-III. 



4.2 Strong phase in D° — > Kit 

The joint decay of D° into K~it + and of D° into a CP eigenstate /„ can be described as 

T Kn , fv ^T[(K-7r + )(f v )} ^A 2 A) v \l + W e- iS \ 2 ^A 2 A 2 fri {l + 2 WC os5), (38) 

where A = | 7r + |?-^|Z) ) | and Af v = \{f rj \'H\D )\ are the real-valued decay amplitudes, r\ = ±1 
is CP eigenvalue of the eigenstate f v , re~ lS is defined in Eq. ([27|) and we have taken f v to be a PP 
or VP CP-eigenstate, without any non trivial phase-space dependence. We also have neglected the 
subdominant r 2 term in Eq. (|38|) . The following asymmetry can be used to determine 5 |10j 

A = l Kw -' U ~ l Kn J - , (39) 

1 KTT-J+ + 1 K-K-J- 
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where ^Kir;f± is defined in Eq. (|38p . which is the rate for the -0(3770) — > D°D° configuration to decay 
into flavor eigenstates and a CP-eigenstates f±. Eq. (|38p implies a small asymmetry, A = 2rcos5. 
In such a case, the error A A is approximately l/\/Nx- 7T +, where Nk- w + is the total number of 
events tagged with CP-even and CP-odd eigenstates, leading to: 

A(cos5)«-^^=. (40) 

The expected number N K - n + of CP-tagged events depends on the total number of D°D° pairs 
N(D°D°), the branching ratio to the CP-eigenstate /„ and the tagging efficiency. Considering all 
decay modes listed in tables [U [2] and [3] we find 

A{cos5)^— r 300 (41) 
^/N(D°D°) 

At BES-III, about 72 x 10 6 D°D° pairs can be collected with four year running |20p29| . which implies 
an accuracy of about 0.03 for cos<5, when considering both K~tt + and K + tt~ final states. 

As in the previous section a more realistic analysis requires a likelihood fit to the full angular 
dependence of the VV modes, which in turn provides independent information on | sin5\ as explained 
above. On the other hand the imperfect reconstruction of the VV events as pure CP-eigenstates will 
presumably introduce sizable systematics in this discussion. 

At a Super r-charm factory [30, 31J with a 2 ab _1 data set, we can expect a factor of ten im- 
provement, but again the precise impact of the modeling of the vector resonances requires more 
studies. 



5 Conclusion 

The charm quark offers interesting opportunities to cross-check the mechanism of CP violation pre- 
cisely tested in the strange and beauty sectors. The start of BES-III will allow for extensive mea- 
surements of charm properties. Among the various tests that can be considered, one may think of 
exploiting the quantum correlations in the DD pairs produced at ^(3770) resonance. In this paper, 
we exploit these correlations in -0(3770) — > D°D° — > (V1V2XV3V4) in connection with CP violation, 
and -0(3770) — > D°D° — > (ViV2)(Kir) for CKM angle 7 measurements, where all VV pairs are 
reconstructed as CP-eigenstates. 

In the case of -0(3770) — > D°D° — > (V1V2XV3V4), the existence of correlations hinders some 
helicity configurations for the outgoing vector mesons in the absence of CP violation. This is mirrored 
by the angular distribution of the differential decay width, out of which CP-violating observables can 
be constructed. Such observables should be interesting to isolate significant New Physics effects in 
the charm sector. Assuming that there would be no CP-violating signal events observed in D°D° 
coherent decays with 20 fb _1 data taken at 0(3770) peak at BES-III, we estimated an order of 
magnitude of the corresponding upper limit on CP-violating parameters, in particular for the channel 
(D°D°)c=-i — > p°p ,K*°p ) ^ . Since the obtained bounds do not follow from a full angular fit 

and do not include the systematics corresponding to the separation of the wanted vector resonances 
from the background, further studies are needed. 

CP-tagged D — > Ktt decays give access to the strong phase difference 8 between Cabibbo-favored 
and doubly-Cabibbo suppressed decays, and thus improve the uncertainty on the 7 measurement of 
the unitary triangle from — > D /DK^ decays. At BES-III, with 20 fb _1 data at -0(3770) peak, we 
estimate the error of cos d to be of a few percents, corresponding to an error on the 5 of a few degrees. 
We expect this estimate can be improved by taking into account the dependence of the full angular 
decay width to the sine of the strong phase. At the Super T-charm factory, the expected statistical 
error on 5 could then fall below one degree. On the other hand a further study of experimental 
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systematics related to the background identification is required since they will presumably dominate 
over the uncertainty quoted here. 

Since our numerical estimates are quite promising we hope that the potential of such coherent 
-D-decays into vector mesons at charm factories will be assessed more precisely in the future. 
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